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ABSTRACT 

We investigate the relationship between HI, H 2 , and the star formation rate (SFR) using azimuthally 
averaged data for seven CO-bright spiral galaxies. Contrary to some earlier studies based on global 
fluxes, we find that Esfr exhibits a much stronger correlation with Eh 2 than with Shi, as Shi saturates 
at a value of ~ 10 M 0 pc -2 or even declines for large Eg FR . Hence the good correlation between SgpR 
and the total (H I+H 2 ) gas surface density S gas is driven by the molecular component in these galaxies, 
with the observed relation between Esfr and Eh 2 following a Schmidt-type law of index n mo i~0.8 if 
a uniform extinction correction is applied or n mo i«1.4 for a radially varying correction dependent on 
gas density. The corresponding Schmidt law indices for S gas vs. SgFR are 1.1 and 1.7 for the two 
extinction models, in rough agreement with previous studies of the disk-averaged star formation law. 
An alternative to the Schmidt law, in which the gas depletion timescale is proportional to the orbital 
timescale, is also consistent with the data if radially varying extinction corrections are applied. We find 
no clear evidence for a link between the gravitational instability parameter for the gas disk ( Q g ) and 
the SFR, and suggest that Q g be considered a measure of the gas fraction. This implies that for a 
state of marginal gravitational stability to exist in galaxies with low gas fractions, it must be enforced 
by the stellar component. In regions where we have both H I and CO measurements, the ratio of H I 
to H 2 surface density scales with radius as roughly i? 15 , and we suggest that the balance between H I 
and H 2 is determined primarily by the midplane interstellar pressure. These results favor a “law” of 
star formation in quiescent disks in which the ambient pressure and metallicity control the formation of 
molecular clouds from H I, with star formation then occurring at a roughly constant rate per unit H 2 
mass. 

Subject headings: galaxies: evolution — galaxies: ISM — stars: formation 


1. INTRODUCTION 

Understanding the factors which control the star forma¬ 
tion rate (SFR) is of fundamental importance to studies 
of the interstellar medium (ISM) and galaxy evolution. 
The most important factor, of course, is the availability of 
cold neutral gas, especially molecular gas. Over the past 
several decades, observations of the CO(l-O) line have es¬ 
tablished a clear link between star formation in the Galaxy 
and molecular gas (e.g., Loren et al. 1973), and even in en¬ 
vironments where star formation is observed without CO 
emission, such as dwarf irregular galaxies, a plausible case 
can be made that molecular clouds are weak CO emit¬ 
ters due to low metal abundances (Wilson 1995; Taylor 
et al. 1998). On the other hand, Kennicutt (1989) has 
shown that the disk-averaged Ssfr is much better cor¬ 
related with Shi than with Eh 2 inferred from CO obser¬ 
vations, a surprising result given the very different radial 
profiles of H I and Ha emission in galaxies. Similar conclu¬ 
sions have been reached by other studies based on global 
averages (Deharveng et al. 1994; Boselli 1994; Casoli et al. 
1996). Thus, it is not obvious whether the H I or H 2 (as 
traced by CO) content will be a better predictor of the 
SFR on large scales. Kennicutt (1989) argues that it is 
the total (H I+H 2 ) gas surface density which shows the 
best correlation with the SFR per unit area. 

A quantitative prescription for the SFR has been pro¬ 
vided by Kennicutt (1998a, hereafter K98), who finds that 
the total gas surface density, averaged over the optical 


disk, is related to the SFR surface density by a Schmidt 
(1959) law, 

Ssfr oc (E gas ) n . (1) 

where n ~ 1.4. Kennicutt’s formulation of the star forma¬ 
tion law has been widely applied in “semi-empirical” mod¬ 
els of galaxy evolution (e.g., Boissier & Prantzos 1999; Tan 
et al. 1999; van den Bosch & Dalcanton 2000). While the 
Schmidt law appears to be valid at high gas surface den¬ 
sities, Kennicutt (1989) and Martin & Kennicutt (2001) 
have presented evidence for a threshold density for star 
formation based on the gravitational instability parameter 
Q. Below the threshold density, which observationally is 
of order 5-10 M 0 pc~ 2 , the Schmidt law breaks down and 
star formation is strongly suppressed. The general appli¬ 
cability of a Q threshold, however, remains controversial 
(Thornley & Wilson 1995; Hunter et al. 1998; Ferguson 
et al. 1998). 

As noted by K98, an alternative formulation of the star 
formation law, where the star formation timescale is pro¬ 
portional to the orbital timescale, 

Usfr oc U ga sU , (2) 

is also consistent with the disk-averaged data. Such a law 
is predicted by models in which spiral arms play an impor¬ 
tant role in triggering star formation (Wyse & Silk 1989), 
or in which star formation is self-regulated to yield a con¬ 
stant value of Q (Larson 1988; Silk 1997). Distinguishing 
between this law and a Schmidt law requires spatially re¬ 
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solved measurements of gas densities and orbital timescales 
(i.e., rotation curves). 

A third prescription for the SFR, arising from CO stud¬ 
ies by J. Young and collaborators (Devereux & Young 
1991; Young et al. 1996; Rownd & Young 1999), states that 
the SFR is roughly proportional to the mass of molecular 
gas, so that their ratio (commonly referred to as the star 
formation efficiency, SFE) is roughly constant. In partic¬ 
ular, Rownd & Young (1999) find that the SFE within the 
disks of spiral galaxies exhibits no strong radial gradients, 
under the assumption of a constant CO-to-H 2 conversion 
factor. An SFE that is independent of galactocentric ra¬ 
dius is also consistent with studies of molecular clouds in 
our Galaxy (Mead et al. 1990; Wouterloot et al. 1988). If 
star formation in disk galaxies can be generically decom¬ 
posed into two processes, the formation of giant molecular 
clouds (GMCs) from H I clouds and the formation of stars 
within GMCs (e.g., Tacconi & Young 1986), then the lat¬ 
ter process occurs at a rate that is largely independent of 
location within a galaxy, and may well be universal in all 
disk galaxies. 

Previous studies of the star formation law in external 
galaxies have tended to concentrate on large samples ob¬ 
served at the low resolution provided by single-dish radio 
telescopes. Although such samples are valuable for char¬ 
acterizing the global properties of galaxies, detailed study 
of individual objects at high spatial resolution is essential 
for understanding the physical basis of the star formation 
law and reconciling the various empirical descriptions. In 
this paper we undertake such an investigation, comparing 
recent CO data from the BIMA Survey of Nearby Galax¬ 
ies (BIMA SONG) with previously published H I and Ha 
data. Our high-resolution CO maps include single-dish 
data, which is essential when accurate flux measurements 
are required for extended sources. We employ azimuthally 


averaged radial profiles in order to improve the signal-to- 
noise ratio and average over time-dependent effects. 

The organization of this paper is as follows. In §2 we de¬ 
scribe our galaxy sample, observing parameters, and data 
reduction procedures. §3 describes how we produced ra¬ 
dial profiles from the intensity images and our corrections 
for Ha extinction. §4.1 presents our results for the relation 
between Esfr and E gas , as well as the relations between 
SgFR and Shi and Eh 2 individually. In §4.2 we compare 
E gas with the threshold surface density discussed by Ken- 
nicutt (1989), and in §4.3 we investigate the relationship 
between the radial profiles of atomic and molecular gas. 
In §5 we show how our results can be interpreted in terms 
of a nearly constant star formation efficiency for molecular 
gas, with molecular cloud formation in turn linked to the 
average ISM pressure. Finally, in §6 we summarize our re¬ 
sults. In subsequent papers of this series, we will analyze 
the CO and H I kinematics of the galaxies in our sam¬ 
ple for evidence of radial inflows (Wong, Blitz, & Bosnia 
2002, Paper II) and evaluate whether the observed gas de¬ 
pletion times and oxygen abundances are consistent with 
closed-box evolution (Wong & Blitz 2002, Paper III). 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Sample Description 

Most of the new observations presented in this paper 
were taken as part of the BIMA Survey of Nearby Galax¬ 
ies (BIMA SONG) (Regan et al. 2001), an interferometric 
CO survey of 44 spiral galaxies selected primarily by op¬ 
tical luminosity (absolute B magnitude Mb < 11). Six 
of these galaxies, along with NGC 4501, were chosen for 
this study. All exhibit strong CO emission at a range of 
radii and have available H I maps from the Very Large Ar¬ 
ray (VLA) at the resolution of the C configuration (~20") 


Table 1 

Properties of the Sample Galaxies 


Name 

RA (2000) 
hh:mm:ss 

DEC 

dd:mm:ss 

Morph. 

Type 

(mag) 

R 25 (b) 

(arcsec) 

(km/s) 

Dist. 

(Mpc) 

NGC 4321 

12:22:54.9 

15:49:21 

SAB(s)bc 

9.98 

220 

1571 

16^) 

NGC 4414 

12:26:27.1 

31:13:24 

SA(rs)c? 

10.62 

110 

716 

ig(e) 

NGC 4501 

12:31:59.2 

14:25:14 

SA(rs)b 

9.86 

210 

2281 

160 

NGC 4736 

12:50:53.1 

41:07:14 

(R)SA(r)ab 

8.75 

340 

308 

4.20 

NGC 5033 

13:13:27.5 

36:35:38 

SA(s)c 

10.21 

320 

875 

17.5^ 

NGC 5055 

13:15:49.3 

42:01:45 

SA(rs)bc 

9.03 

380 

504 

10(5) 

NGC 5457 

14:03:12.5 

54:20:55 

SAB(rs)cd 

8.21 

870 

241 

7 . 4 O 


“Total blue magnitude, corrected for Galactic and internal absorption and redshift, from RC3 (de Vaucouleurs et al. 1991). 
^Semimajor axis at 25.0 mag arcsec -2 , from RC3. 

“Heliocentric velocity, from NED. 

d Cepheid distance to NGC 4321 from Ferrarese et al. (1996) adopted for Virgo cluster. 

“Cepheid distance from Turner et al. (1998). 

^Redshift distance using Hq = 75 km s -1 Mpc -1 . 

^Distance from Mark III data (Willick et al. 1997) using Ho=75. 

^Cepheid distance from Kelson et al. (1996). 
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or better. The insistence on high-resolution H I data is 
essential for the CO and H I data to be used in a com¬ 
plementary way; although H I maps at even higher (~6") 
resolution would be ideal, few such maps are available be¬ 
cause of the large investment of observing time needed to 
obtain good sensitivity at those resolutions. The galax¬ 
ies span a range of morphological type, from Sab to Scd. 
While none are involved in strong interactions with other 
galaxies, NGC 4321 (M100) and NGC 4501 are members 
of the Virgo Cluster and presumably interacting with the 
intracluster medium. The basic properties of the sample 
are summarized in Table 1. The tabulated coordinates are 
for the center of the galaxy as given by the NASA Extra- 
galactic Database (NED), except for NGC 5055 where we 
have adopted the peak of the astrometrized AT-band image 
of Thornley & Mundy (1997a). 
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Fig. 1.— Histogram of H 2 masses normalized by B luminosity, 
using data from Young et al. (1995) and a constant X-factor, for our 
sample ( shaded ) and 40 of the 44 galaxies in the optically selected 
BIMA SONG sample. The remaining 4 galaxies were not observed 
by Young et al. (1995) and probably have below-average CO fluxes. 
For NGC 5457 (second shaded box from left) we have used the CO 
flux measured by Kenney et al. (1991), which is probably an under¬ 
estimate since their map is incomplete. The vertical dashed line is 
for the Milky Way. 



Based on single-dish measurements from the Five Col¬ 
lege Radio Astronomy Observatory (FCRAO) Extragalac- 
tic CO Survey (Young et al. 1995), the median CO flux 
for the galaxies in our sample is 2760 Jy km s^ 1 , com¬ 
pared to a median of 1640 Jy km s -1 for the entire SONG 
sample. To determine if the difference results from our 
choice of more optically luminous galaxies, Figure 1 shows 
a histogram of the ratio of H 2 mass to R-band luminosity, 
as tabulated by Young et al. (1995), for the galaxies in 
our sample (shaded) as well as for the entire SONG sam¬ 
ple. We have scaled the CO fluxes derived by Young et al. 
(1995) to H 2 masses using a conversion factor appropriate 


for the Galaxy: 

% = 2 x 10 20 


Ico 


cm 


( 3 ) 


K km s 1 

(Strong & Mattox 1996; Dame et al. 2001). The median 
for the entire SONG sample (—1.17 dex in solar units) is 
somewhat less than the value for the Milky Way (—1.1 
dex, based on Dame 1993 and van der Kruit 1990), which 
in turn is less than the median for our sample (—1.0 dex). 
Note that a difference of 0.17 dex corresponds to a factor 
of 1.5. Despite considerable uncertainty in the FCRAO 
fluxes, since they are estimated from major axis scans 
rather than complete two-dimensional maps, we conclude 
that the galaxies in our sample are somewhat overluminous 
in CO (relative to their optical luminosities) compared to 
bright galaxies in general. This is not surprising, given 
that we explicitly required that CO be observed over a 
large portion of the galaxy disk. Thus, one should bear 
in mind that properties of the ISM in the galaxies studied 
here may not be representative of the field galaxy popula¬ 
tion, or even the population of optically luminous galaxies. 


2.2. BIMA CO Observations 

Observations with the BIMA 2 interferometer at Hat 
Creek, California were conducted primarily in the C ar¬ 
ray configuration, which gives a synthesized beam of ~6" 
FWHM, with a few additional observations (“tracks”) 
taken in the D array configuration (15" FWHM). The syn¬ 
thesized beam gives the resolution of the resulting maps, 
and is determined by the maximum spacing of the anten¬ 
nas: in C-array the spacings range from 7 to 90 m (2.7-35 
kA), while in D-array the spacings range from 7 to 30 m 
(2.7-11 kA). Several of the C-array observations of NGC 
4501, 4736, and 5033 were made in 1996 and 1997 when 
the BIMA array consisted of nine antennas; all other ob¬ 
servations were made in 1998 and 1999 with ten antennas. 

In our standard frequency setup, the receivers were 
tuned to the CO (J = 1 —> 0) transition at 115.2712 GHz 
(A=2.6 mm), with the correlator configured to have 4 in¬ 
dependently positioned spectral windows, each with 100 
MHz bandwidth and 64 channels. The spectral windows 
were placed side-by-side in frequency with an overlap of 
«10 MHz (6 channels). The resulting velocity range cov¬ 
ered was «950 km s -1 at a resolution of 4.06 km s -1 . Data 
from the lower sideband of the first local oscillator, at a 
frequency of 112 GHz, were also recorded, but no contin¬ 
uum emission was detected in any of the sample galaxies. 

At A=2.6 mm the BIMA antennas have a primary beam 
FWHM of 100", which defines the effective field of view. 
During most observing sessions a 7-point hexagonal grid 
centered on the galaxy nucleus was observed, with the 
points separated by 50" (for NGC 4501 and 4736) or 44" 
(for the remaining galaxies). In order to maintain roughly 
equal sampling in the visibility plane for each pointing cen¬ 
ter, the entire grid was observed in 7-9 minutes and then 
repeated. For a few tracks on NGC 4414, 4736, and 5033, 
only the central pointing was observed. In cases where the 
same galaxy was observed as both a single pointing and a 
mosaic, the tracks were imaged together as a mosaic—with 
the exception of NGC 4414, where only the central point¬ 
ing was used for the high-resolution map. The resulting 
values of i? max for each datacube, defined as the radius at 
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Table 2 

Properties of the BIMA CO maps 


Galaxy 

Weighting 

Resolution 
(" x " x km s" 1 ) 

PAW 

(°) 

(mjy/bm) 

R V. 

1 l max 

(") 

NGC 4321 

natural 

8.85 

X 

5.62 x 10 

4.7 

49 

94 


20" taper 

20 

X 

20 x 10 

0.0 

97 


NGC 4414 

robust 

6.53 

X 

4.88 x 10 

8.2 

42 

50 


13"taper 

17.5 

X 

15.3 x 10 

72.1 

64 

94 

NGC 4501 

natural 

7.89 

X 

6.23 x 20 

-80.9 

56 

100 


22" taper 

23.2 

X 

22.3 x 10 

24.6 

174 


NGC 4736 

robust 

6.86 

X 

5.02 x 10 

62.5 

69 

100 


13"taper 

15 

X 

15 x 10 

0.0 

83 


NGC 5033 

natural 

7.52 

X 

6.49 x 20 

82.8 

52 

94 


18"taper 

19.9 

X 

17.0 x 20 

89.1 

108 


NGC 5055 

robust 

5.69 

X 

5.30 x 10 

13.9 

83 

94 


12"taper 

12.9 

X 

12.8 x 10 

-83.5 

107 


NGC 5457 

natural 

6.88 

X 

6.36 x 10 

72.5 

43 

94 


15" taper 

15 

X 

15 x 10 

0.0 

80 



“Position angle of beam measured from N towards E. 

6 RMS noise in a single emission-free channel measured across the inner 2' X 2', in mjy bm _1 . 
“Radius at which half-power point of the outermost pointing is reached. 


which one reaches the half-power point of the outermost 
pointing, are given in Table 2. 

2.2.1. Assessment of Data Quality 

The two principal factors which can affect the data qual¬ 
ity from a millimeter-wave interferometer are the tropo¬ 
spheric water vapor column, which increases the atmo¬ 
spheric opacity, and fluctuations in that column, which 
lead to variations in the refractive index of the atmosphere 
(“seeing”) and hence variations in the phase of the signal 
measured at the telescope. Changes in atmospheric opac¬ 
ity are reflected in the “effective” system temperature T* , 
which ranged from 300-600 K during our observations, de¬ 
pending on the source elevation and weather conditions. 
Atmospheric phase noise can be quantified by a frequency- 
independent quantity a pa th, representing the RMS fluctu¬ 
ation in the difference in effective path length through the 
atmosphere for two parallel sight lines separated by a given 
distance (e.g., Wright 1996). Since mid-1998, (J pa th has 
been measured every 10 minutes at Hat Creek by a ded¬ 
icated phase monitor, which records signals from a tele¬ 
vision satellite using a pair of commercial satellite dishes 
separated by 100 m (Lay 1999). These measurements yield 
a fiducial cr pa th,o , which was typically 200-400 /im for the 
C-array observations and 800-1200 /im for the D-array ob¬ 
servations (the latter occurring mostly in summer). Note 
that with the RMS phase given by 4> rm s = 2-7T cr pa th/\ the 
visibility amplitude after time-averaging is given by 

(V) = V 0 e~ ( ' 4 ’ rms)2/ ‘ 2 . (4) 

Hence a value of cr pa th,o = 560 /im corresponds to a 50% 
decorrelation (i.e. loss of visibility amplitude) at A = 3 
mm for a 100 m baseline (Akeson 1998), roughly the 


longest baseline in the C-array. Although the measured 
(jpath ,o was typically worse than this in D-array, the ac¬ 
tual phase fluctuations that affect the short baselines of 
D-array (b <C 100 m) are much smaller since a pa th oc b 5 ' e 
(Wright 1996). 

2.2.2. Data Calibration 

The data were calibrated and reduced following the pro¬ 
cedures described in Regan et al. (2001). We give only a 
summary here. An online linelength calibration was ap¬ 
plied to correct for variations in the electrical path length 
from the antenna to the correlator. Remaining drifts in 
the antenna phase gains were then corrected by observ¬ 
ing a nearby (within ~15°) quasar every half hour and 
fitting a low-order (quadratic or cubic) polynomial to the 
antenna-based self-calibration solutions. Mars was used 
for amplitude calibration except in 1998 April, when the 
quasar 3C273 was used. For all tracks both amplitude and 
phase gains were assumed to be frequency-independent; 
this approximation was found to be adequate in several 
cases by observing 3C273 using the same correlator setup 
as for the source observations. 

For tracks during which a pa t.h was measured, an addi¬ 
tional calibration step was performed in order to compen¬ 
sate for the effects of phase decorrelation on the maps. 
In this procedure, implemented in a MIRIAD task called 
UVDECOR written by M. Regan, the visibility amplitudes 
for all datasets are scaled up by a baseline-dependent cor¬ 
rection factor, typically ~1 1.2, based on the expected am¬ 
plitude loss due to phase decorrelation (see Regan et al. 
2001). We found that applying UVDECOR led to a ~10% 
increase in the fluxes of the final maps. Overall, we esti¬ 
mate that the flux scale of the resulting BIMA maps is 
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accurate to within ~20%. 

Following the initial calibration, the visibilities were 
subjected to three iterations of phase-only self-calibration, 
which in all cases led to convergence. Channel maps were 
then produced for each galaxy by assigning weights to the 
calibrated visibilities as described below, transforming to 
the image plane after binning to 10 or 20 km s -1 chan¬ 
nels (to improve signal-to-noise), and deconvolving using 
a variant of the CLEAN algorithm developed by Steer 
et al. (1984) to improve the deconvolution of spatially ex¬ 
tended emission. In general, assigning greater weight to 
short-baseline data results in lower resolution but higher 
signal-to-noise for extended sources, as well as reduced sys¬ 
tematic errors because short baselines are less subject to 
atmospheric phase noise. Thus, for all galaxies we gener¬ 
ated a low-resolution map by applying a Gaussian taper 
in the visibility plane, chosen to match the resolution of 
the galaxy’s H I cube. In addition, a high-resolution map 
was generated using either natural or robust weighting. 
Natural weighting provides the best point-source sensitiv¬ 
ity and was used for sources with poorer signal-to-noise 
(NGC 4321, 4501, 5033, and 5457). Robust weighting, 
which provides poorer sensitivity but somewhat higher res¬ 
olution, was used for NGC 4414, 4736, and 5055. Table 2 
summarizes the properties of the final deconvolved maps. 

2.3. Kitt Peak CO Observations 

Single-dish CO observations were performed at the 
NRAO 3 12m telescope on Kitt Peak, Arizona. We refer 
to these data hereafter as the “KP data.” At A=2.6 mm 
the half-power beamwidth (HPBW) of the telescope is 
55". All galaxies with the exception of NGC 4501 were 
observed between 1998 April and 2000 January using the 
On-The-Fly (OTF) technique, in which the telescope is 
scanned rapidly across the sky and spectra are written 
out at a very high rate (every 0.1 s). The relative point¬ 
ing of the telescope during the map is also continuously 
measured, although the absolute pointing must still be 
checked and adjusted at regular intervals (once every 1 
2 hr) by observing a bright continuum source (typically 
Mars). Each map covered a field roughly 6' x 6 ' in size and 
was completed in <15 minutes; roughly 20-40 maps were 
obtained per galaxy and later co-added during data reduc¬ 
tion. NGC 4501 was observed in 1997 June in position- 
switched (hereafter PS) mode, where spectra are taken at 
specified grid points spaced by half a beamwidth or less. 
The integration time was typically 4-8 minutes per point¬ 
ing, with a total of 27 pointings observed covering a region 
roughly 2! x 2 ' in size. 

The telescope was configured to observe orthogonal 
polarizations using two 256-channel filterbanks (one for 
each polarization) at a spectral resolution of 2 MHz (5.2 
km s _1 ). Calibration was performed using the chopper 
wheel method every 5-15 minutes; typical values of T* ys 
were 300-400 K. Spectra were output on the T* R scale, 
which represents the source antenna temperature cor¬ 
rected for ohmic losses, atmospheric absorption, spillover, 
and scattering (Kutner & Ulich 1981). Background emis¬ 
sion was measured by observing an off-source position 10' 
from the pointing center (for the PS observations) or by 


averaging spectra at the edges of the OTF maps; the off- 
source spectrum nearest in time was then used for back¬ 
ground subtraction. 

Data reduction was performed in the COMB software 
package for the PS spectra and the AIPS and MIRIAD 
packages for the OTF spectra. A linear baseline was fit¬ 
ted to the signal-free channels of each spectrum and sub¬ 
tracted. FITS datacubes were then generated using a cell 
size of 13'.'5 for the PS map and 18" for the OTF maps. 
For each galaxy, the RMS noise in the final maps ranged 
from 14-28 mK in a 2 MHz channel. The maps were con¬ 
verted into Jy units using a standard telescope gain of 33 
Jy (KpS])" 1 . 

Prior to combination with the BIMA maps, OTF maps 
from different dates were cross-correlated against each 
other to correct for pointing offsets, as described in Regan 
et al. (2001). The derived offsets were generally around 
10" but ranged as high as 27" (half a beam width) in one 
case. After this internal registration of the OTF maps, 
a final cross-correlation with the BIMA map (smoothed 
to 55") was performed, under the assumption that the 
flux resolved out by BIMA did not significantly affect the 
centroid of the CO emission. The resulting corrections 
amounted to 5"-ll". For NGC 5457, which has a low in¬ 
clination and a large angular diameter, this assumption 
was found to be inappropriate due to the large amount 
of flux missing from the BIMA map, and the final cross¬ 
correlation was not performed. 

2.4. CO Data Combination 

The relative calibration between two millimeter tele¬ 
scopes is always somewhat uncertain, so when combining 
single-dish and array data one typically attempts to solve 
for a relative calibration factor, using data from both tele¬ 
scopes that have been filtered to the same range of spatial 
frequencies. The factor needed to scale the KP data to 
match the BIMA data was calculated for each galaxy us¬ 
ing the MIRIAD task IMMERGE and ranged from 1 to 
1.3. In the case of NGC 4501, a reliable calibration factor 
could not be determined since the field of view of the KP 
observations was not large enough, so a factor of 1 was 
adopted. We note that it is quite possible that the BIMA 
data should be scaled down rather than the KP data scaled 
up; this is especially likely if the effect of decorrelation 
has been overestimated when applying UVDECOR. For 
simplicity, we assume here that the BIMA data have the 
“correct” flux scale and emphasize that the CO fluxes we 
measure may be slightly overestimated. Given the uncer¬ 
tainties in converting CO brightness to H 2 column density, 
an overall scaling factor of <1.3 is probably unimportant. 

We examined two methods of combining the BIMA and 
KP data, one which generates a weighted sum of the two 
cubes before deconvolution (the “dirty-map” combination 
method) and another which folds in the KP data after 
the BIMA map has been deconvolved, by partially replac¬ 
ing the inner part of the Fourier plane with the KP data. 
The two methods gave similar results to within the noise. 
We adopted the latter method, as implemented in the IM¬ 
MERGE task. The resulting datacubes are sensitive to 
structures on all spatial scales larger than ~6". A more 


3 The National Radio Astronomy Observatory is a facility of the National Science Foundation, operated under cooperative agreement by 
Associated Universities, Inc. 
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Table 3 

Properties of the VLA HI Maps 


Galaxy 

Resolution 
(" x " x km s -1 ) 

PA+) 

(°) 

(rnjy/bm) 

Reference 

NGC 4321 

20 x 20 x 20.6 

0.0 

0.75 

Knapen et al. (1993) 

NGC 4414 

17.5 x 15.3 x 10.3 

72.1 

0.4 

Thornley & Mundy (1997b) 

NGC 4501 

23.2 x 22.3 x 10.3 

24.6 

0.6 

A. Bosma (Paper II) 

NGC 4736 

15 x 15 x 5.15 

0.0 

1.8 

Braun (1995) 

NGC 5033 

19.9 x 17.0 x 20.6 

89.1 

0.3 

Thean et al. (1997) 

NGC 5055 

12.9 x 12.8 x 10.3 

-83.5 

0.2 

Thornley & Mundy (1997a) 

NGC 5457 

15 x 15 x 5.15 

0.0 

1.8 

Braun (1995) 


“Position angle of beam measured from N towards E. 

6 RMS noise in a single emission-free channel measured across the inner 4' X 4', in mjy bm _1 . 


detailed discussion of issues surrounding the combination 
of single-dish and interferometer data can be found in Ap¬ 
pendix B of Wong (2000). 

2.5. HI Data 

The H I data used in this study were obtained at the 
NRAO Very Large Array (VLA) by other researchers, who 
kindly provided us with reduced clatacubes in FITS for¬ 
mat. The properties of these cubes and references to the 
original papers in which they were presented are given in 
Table 3. (Details of the NGC 4501 observations will be 
presented in Paper II.) At A=21 cm the HPBW of the 
VLA is 32'. All galaxies were observed in the C configu¬ 
ration, which gives a synthesized beam of ~13 r/ for uni¬ 
form weighting and ^20” for natural weighting; three of 
the galaxies (NGC 4736, 5055, and 5457) were also ob¬ 
served at higher resolution in the B configuration. Data 
obtained in the D configuration are included for all galax¬ 
ies except NGC 4501 and 5033, so that the resulting maps 
are sensitive to even very extended (~10 / ) emission. The 
datacubes were imported into MIRIAD and transformed 


to the spatial and velocity coordinate frames of the BIMA 
cubes. 

2.6. Ha Images 

Continuum-subtracted Ha+[N II] CCD images were ob¬ 
tained from a variety of sources, as listed in Table 4. All 
images were originally taken at the Kitt Peak National Ob¬ 
servatory (KPNO) 0.9-m telescope, except for NGC 5033 
(Palomar 60-inch) and 5457 (KPNO Burrell Schmidt tele¬ 
scope). For NGC 4321 and 5055, the images were obtained 
in 1999 April under photometric conditions, with the stan¬ 
dard stars Feige 34 and PG 0939+262 used for flux cali¬ 
bration (Massey et al. 1988). The other images have not 
been flux calibrated, so the flux scale was determined by 
comparison with published Ha fluxes in the literature, as 
stated in the table. Note that the fluxes of Kennicutt & 
Kent (1983) have been scaled up by a correction factor of 
1.16 (K98). The uncertainty in our bootstrap calibration 
is at least as large as the uncertainty in the original pub¬ 
lished fluxes (generally ~25%). However, this will only 
affect the overall flux level, which is subject to significant 


Table 4 

Sources for Ha Images 


Galaxy 

Image 

Size 

- log /hq 

1°> Credit 

Photometry 

Reference 

NGC 4321 

IT 

11.28 

M. Regan, K. Sheth 

Massey et al. (1988) 

NGC 4414 

IT 

11.38 

Thornley & Mundy (1997b) 

Thornley & Mundy (1997b) 

NGC 4501 

6 ' 

11.14 

Martin & Kennicutt (2001) 

Young et al. (1996) 

NGC 4736 

6 ' 

10.80 

Martin & Kennicutt (2001) 

Kennicutt & Kent (1983) 

NGC 5033 

12 ' 

11.29 

M. Regan 

Kennicutt & Kent (1983) 

NGC 5055 

IT 

10.74 

M. Regan, K. Sheth 

Massey et al. (1988) 

NGC 5457 

20 ' 

10.18 

van Zee et al. (1998) 

Kennicutt (1988) 


a fna i s the total flux in the image, given in erg cm 2 s 1 . 
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uncertainties anyway when converted into a star formation 
rate (§3.2.3). 

Since most of the images contained artifacts due to bad 
pixels, cosmic ray events, and foreground stars that were 
not cleanly subtracted (e.g., due to changes in seeing), 
these were masked out interactively using the IRAF task 
IMEDIT. Artifacts isolated from emission were replaced 
with zero values, while those superposed on diffuse emis¬ 
sion were replaced by a background substitution algorithm 
that fits a smooth surface to an annular region of the image 
centered on the artifact. Remaining outliers were removed 
by clipping negative values below — 5er and positive val¬ 
ues above a limit determined by visually identifying the 
brightest H II region. Astrometry was performed by com¬ 
paring the locations of stars in the continuum image with 
a DSS image using the KARMA package (Gooch 1995) 
and then transferring the plate solution to the Ha image; 
in some cases there were compact Ha sources visible in 
the DSS image and the astrometry could be determined 
by direct comparison to the Ha image. Astrometric errors 
should be less than 2". 

3. DATA ANALYSIS 
3.1. Intensity Images 

Integrated intensity (“moment-0”) images were pro¬ 
duced from the CO and H I datacubes using a masking 
technique designed to exclude noise. The necessity for 
such a procedure arises from the fact that the emission 
at any given pixel is usually confined to a small range of 
velocities, so a direct sum of all channels in the datacube 
would result in very poor signal-to-noise (see discussion in 
Bosnia 1981; van der Kruit & Shostak 1982). Thus, em¬ 
ploying a similar technique to Braun (1995), we defined a 
search region in each channel map by the 3cr smo contour of 
a smoothed version of the map, and included pixels within 
that region with intensities |/| > tr c h in the sum (here <r smo 
and <T c h are the RMS noise per channel in the smoothed 
and original cubes respectively). Tests showed that while 
use of this masking technique does not have a significant 
effect on the total flux in the datacube, it does improve 
the accuracy with which weaker emission features are rep¬ 
resented. 

The resulting intensity images are displayed in Figure 2. 
For each galaxy, the top panels show the CO intensity map 
at full resolution and the corresponding Ha image for com¬ 
parison. The field of view of the BIMA mosaic (where the 
sensitivity drops to half of its peak value) is shown as a 
dotted contour. Most galaxies show a strong correspon¬ 
dence between CO and Ha emission, although the emission 
peaks do not always coincide. In this respect, NGC 4736 
appears unusual since the CO and Ha distributions show 
fundamentally different symmetries, as discussed in Wong 
& Blitz (2000). With the exception of NGC 4736 and 5033, 
the CO emission appears to extend to the edge of the ob¬ 
served field, suggesting that our observations sample only 
the inner part of the molecular disk. (At a mean distance 
of 13 Mpc, our CO observations extend to R ~ 6.3 kpc.) 

The bottom panels for each galaxy show the tapered 
CO intensity map alongside the VLA H I intensity map. 
Both maps have the same spatial resolution (indicated by 
the beam in the lower left corner of each panel). We find 
that four of the galaxies (NGC 4321, 4501, 4736, 5457) 


show both an H I deficit and CO enhancement in their in¬ 
ner regions, suggesting a large-scale conversion of atomic 
to molecular gas. This phenomenon has long been noted 
(e.g., Morris & Lo 1978), and its occurrence appears to be 
unrelated to the Hubble type of the galaxy, since 4736 is 
classified as early-type (Sab) whereas 5457 is classified as 
late-type (Scd). 


3.2. Radial Profiles 


Radial profiles were derived from the intensity images 
by averaging in elliptical annuli spaced by roughly half 
the synthesized beamwidth and centered on the position 
given in Table 1. The annuli were assumed to be circular 
rings viewed at an inclination i and with a line of nodes 
rotated from north to east by some position angle (PA); 
the adopted values of PA and i for each galaxy are given in 
Table 5. As described in Paper II, the appropriate PA for 
each galaxy was determined by kinematic fits to the ve¬ 
locity field, while i was determined from isophotal and/or 
kinematic fits. The mean intensity for each annulus was 
assigned an uncertainty of 


^mean 


J pix 


\J N pix/'ft'be 


( 5 ) 


where cr p i x and IVpi x are the RMS pixel value and the num¬ 
ber of pixels in the annulus, and nbeam is the number of 
pixels in each resolution element (“beam”). The resolu¬ 
tion of the Ha images was assumed to be 2" (Gaussian 
FWHM). Note that <7 mean is a quantity which reflects only 
the dispersion in the values within each annulus, and thus 
reflects a combination of measurement errors and any in¬ 
trinsic non-axisymmetry. 


Table 5. Adopted Orientation Parameters 


Galaxy 

P.A>) 

(deg) 

(deg) 

NGC 4321 

153 

34 

NGC 4414 

159 

55 

NGC 4501 

140 

63 

NGC 4736 

295 

35 

NGC 5033 

353 

68 

NGC 5055 

98 

63 

NGC 5457 

42 

21 


“Position angle of receding side of line of nodes, measured E from 
N. 


3.2.1. CO and HI Profiles 


Since the CO and H I intensity images are generated 
from datacubes by a masking process, sampling effects 
make their noise properties somewhat unusual. Especially 
at large radii where few pixels survive the masking, the 
RMS fluctuation within an annulus (given by cr mea n above) 
can greatly underestimate the error in the profile because 
masked pixels have been assigned a value of zero. We 
can derive a conservative estimate of the true uncertainty 
in the mean for each annulus by considering the case in 
which all channels are included in the moment map (i.e., 
no masking or clipping of any kind): 


^est — 


(Teh An y/nffi 
\J Api x / Rbe a m 


( 6 ) 





NGC 4321 


NGC 4414 

5D 


Fig. 2.— Intensity images for each of the seven galaxies in our sample. Top left: CO at full spatial resolution (robust or natural weighting). 
Contours are kin 2 Jy bm -1 km s _1 , for integers n > 2 and k±=0.75 for NGC 4321, 4414, and 4501, k\ = l for NGC 4736 and 5055, ki=1.5 
for NGC 5033, and /ci=0.6 for NGC 5457. Top right: Ha image with CO contours overlaid. Bottom left: CO at the resolution of the H I 
image. Contour levels are k^n 2 Jy bm -1 km s -1 , for integers n > 2 and &2=2 for NGC 4321, 4414, 4501, and 4736, &2=4 for NGC 5033, 
&2=3 for NGC 5055, and &2=1 for NGC 5457. Bottom right: H I intensity image, with contours of 0.02n 2 Jy bm -1 km s -1 (n > 2). 
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Fig. 3.— Radial profiles of H I, H 2 , and SFR surface density. Top panels: H I and H 2 profiles, corrected for inclination. The H 2 profile 
derived from the BIMA CO data alone is shown as a dashed line. Unfilled plot symbols represent uncertain measurements that may suffer 
from incompleteness (see text). Bottom panels: SFR profile (triangles), corrected for inclination but not for Ha extinction. The short and 
long dashed lines are the resulting SFR for uniform and TVh- dependent extinction corrections respectively. 
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where cr c h is the RMS noise in a channel map of the orig¬ 
inal datacube, n c h is the number of channels, and Av is 
the velocity width of each channel. 

The derived radial gas profiles, at the resolution of 
the H I maps, are presented in the top panels of Fig¬ 
ure 3. Points are separated by 8", oversampling the pro¬ 
files slightly. Unfilled plot symbols are used to denote 
radii where the mean intensity falls below 3<j es t, since these 
measurements are less secure, relying on the ability of the 
mask to select out the “real” emission. These radii are 
also the most likely to suffer from incompleteness, which 
can result if the mask excludes low-level emission that still 
makes a substantial contribution to the sum when inte¬ 
grated over a large annulus. The profiles are given in mass 
surface density units, which are derived as follows. First 
the conversion from surface brightness to column density 
is performed using the CO-to-kU conversion factor given 
by Equation 3, or using the optically thin approximation 
for HI: 


7V H i = 1.82 x 10 18 


K km s 


-l 


cm 


-2 


(7) 


where I 21 is the velocity-integrated H I intensity. Column 
densities were then converted to face-on surface densities, 


Ehi = wh Nm cos* , (8) 

Eh 2 = 2mn 1 Vh 2 cos i , (9) 

excluding helium. The BIMA-only CO profiles, converted 
to the appropriate units, are also shown for comparison. 
Note that the extra flux provided by the KP data has lit¬ 
tle effect on the CO radial profiles, aside from a constant 
offset which only becomes significant at outer radii. 

The assumption of a constant CO-to-H 2 conversion fac¬ 
tor (A'-factor) remains quite controversial, with some re¬ 
searchers suggesting that X is increased by a factor of 10 
or more in the vicinity of the Galactic Center (Sodroski 
et al. 1995; Dahmen et al. 1998). Indeed, the presence 
of central CO condensations in many galaxies (Sakamoto 
et al. 1999; Regan et al. 2001) may be an indication of a 
change in the A'-factor rather than a true concentration 
of molecular gas. On the other hand, two of the seven 
galaxies (NGC 4501 and 5033) show featureless, roughly 
exponential CO profiles. We discuss the possibility of a 
varying A-factor further in §5.1. 


3.2.2. Total Gas Profiles 

Radial profiles of the total gas surface density (E gas ) 
were derived by summing the contributions from H I and 
H 2 and multiplying a factor of 1.36 to include helium. 
Since the CO emission is only imaged within R < 2', and 
our sensitivity to CO emission falls off at H 2 column den¬ 
sities comparable to the peak H I column densities of ~10 
M 0 pc -2 , it was necessary to extrapolate the CO profile 
to larger radii assuming an exponential profile. The scale 
length for the extrapolation was determined from a fit to 
the entire radial profile, except for NGC 4414, where the 
fit was made to the declining part of the CO profile, and 
NGC 5457, where the scale length was determined from 
the radial profile of a large (10' x 10') CO map obtained 
at the NRAO 12m telescope as part of BIMA SONG. To 
take an example, Figure 4 shows the total gas profile for 
NGC 4321 and its individual components (also corrected 
for helium). The profiles are probably most uncertain at 


radii just beyond i? max (refer to Table 2), where the ex¬ 
trapolated H 2 mass is likely to be significant. Even so, we 
expect molecular gas to contribute less than half of the 
gas mass at these radii, and hence the error in S gas in¬ 
troduced by the extrapolation should be <50% as long as 
the CO profile continues to decline beyond i? max , an as¬ 
sumption supported by visual inspection of our somewhat 
larger OTF maps. 


NGC 4321 



Fig. 4.— Total gas profile for NGC 4321 derived from the CO 
and H I profiles, with the CO profile extrapolated as shown. The 
error bar in the upper left corner represents a factor of 2 change in 
any of the curves. 


3.2.3. Ha Profiles 


The Ha emission from spiral galaxies is dominated by 
young {t < 20 Myr), massive (>10 M 0 ) stars, and is thus 
a reasonable tracer of the recent star formation rate (SFR) 
(see the review by Kennicutt 1998b). We expressed the Ha 
brightnesses in units of SFR density using the conversion 
(Kennicutt 1998b): 


SFR 


= 7.9 x 10~ 42 


J Ha 


( 10 ) 


M 0 yr -1 \ergs -1 

This conversion ignores extinction by dust and assumes a 
Salpeter IMF between 0.1 100 M 0 , a constant SFR with 
time, and negligible escape of ionizing radiation from the 
galaxy. Since the Ha images include a contribution from 
[N II], they have been multiplied by 0.67, assuming an 
average ratio of [N II] /Ha ~ 0.5 (Kennicutt 1992). Al¬ 
though the use of a uniform correction is rather simplistic, 
since [N II]/Ha is often enhanced in the nuclear region 
(e.g., Rubin & Ford 1986), it should be roughly applicable 
for most of the H II regions in the disk. For comparison, 
Smith et al. (1991) used long-slit spectroscopy to deduce a 
ratio of Ha/([N II]+Ha) « 0.62 for the star-forming ring 
in NGC 4736. 


The bottom panels of Figure 3 show the radial Ha pro¬ 
files, converted to units of SFR density, derived from im¬ 
ages smoothed to the same resolution as the H I images. 
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These can be directly compared to the gas profiles in the 
top panels. The data points have not been corrected for 
extinction, whereas the dashed lines represent the results 
of applying the extinction corrections described in §3.3. 
In §4.1.1 we also compare the CO and Ha profiles at the 
higher resolution of the CO data, although we have not 
shown the corresponding profiles here. 

3.3. Extinction Corrections for Ha 

While extinction at Ha due to Galactic dust along the 
line of sight is probably unimportant (Ay^ai < 0.13 mag 
for all seven galaxies according to Schlegel et al. 1998), 
extinction internal to the galaxies themselves is likely to 
be significant, given that stars form in regions of dense gas 
that have associated dust. Thus, when deriving quantita¬ 
tive estimates of the SFR from Ha emission, it is important 
to take extinction into account. 

In his study of disk-averaged star formation, K98 ap¬ 
plied a constant extinction correction of 1.1 mag to his 
Ha- data, based on an earlier comparison of integrated Ha 
and radio continuum fluxes (Kennicutt 1983). Although 
obviously a gross simplification, it is difficult to do better. 
The use of an Ha/H [3 ratio to determine the extinction is 
inapplicable for galaxy-integrated spectra, due to the dif¬ 
ficulty of measuring the H f3 flux in the presence of stellar 
absorption lines (Kennicutt 1992). Even with spectropho¬ 
tometry of individual H II regions, one can typically only 
obtain measurements for the brightest sources, and the de¬ 
rived extinctions may not be appropriate for all of the Ha 
emission (Bell & Kennicutt 2001). Hence, the most reli¬ 
able extinction estimates are likely to come from narrow- 
band imaging of NIR recombination lines such as Bry or 
Paa, or from measurements of thermal radio continuum 
fluxes. 

We have considered four techniques for determining ex¬ 
tinction corrections for our sample: 


1. One can apply a global correction following K98, 
but instead using measured radio continuum and Ha 
fluxes for individual galaxies. Niklas et al. (1995) 
published radio fluxes at 10 GHz for a large set of 
galaxies that encompasses our sample, and Niklas 
et al. (1997) have combined these with data at other 
frequencies to separate the thermal and non-thermal 
components (the latter often being dominant even 
at 10 GHz) via the spectral slope. The thermal flux 
Sth can be converted into a predicted Ha flux using 
the formula (Condon 1992): 

F p ,, d (H„) = 8.0x10-'-’ ( g ^)° 1 erg cm'’ 

(n) 

This relation assumes an electron temperature of 
10 4 K, but is not very sensitive to it (-F pre d a T~ 0 5 ). 
In Figure 5(a) we compare the integrated fluxes 
in our Ha images with the radio-predicted fluxes. 
While extinction appears to be negligible for NGC 
5457 (M101), a typical value of 4 ho ~ 1-5 mag is 
found for the other six galaxies. For these galaxies, 
the SFR is therefore underestimated by a factor of 
~4 when using the Ha fluxes alone. 

2. One can derive a radially dependent correction us¬ 
ing the observed radial gas profile and an assumed 
dust-to-gas ratio. We first assume a “foreground 
dust screen model,” in which half of the observed 
gas column density is located in a uniform absorb¬ 
ing slab between the Ha emission and us. We adopt 
standard Galactic values for the gas-to-dust ratio, 
Nyi/Av = 2 x 10 21 cm -2 mag -1 , and the extinc¬ 
tion law, Ar/Av = 0.75 (Bohlin et al. 1978; Rieke 
& Lebofsky 1985). The inferred extinctions reach 
peak values of 4 14 mag near the galaxy center, al¬ 
though even these represent azimuthally averaged 
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continuum data. Dashed lines correspond to extinctions of 0 and 1.5 magnitudes, (b) Comparison of the disk-averaged extinction at Ha 
derived from comparing the radio and Ha fluxes with the extinction predicted by a “hybrid” dust model, where dust associated with H 2 is 
well-mixed but dust associated with H I constitutes a foreground screen. The model assumes a Galactic dust-to-gas ratio and extinction law. 
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values at ~15" resolution. After correction for ex¬ 
tinction, the global SFRs increase by factors ranging 
from 1.3 (for NGC 5457) to 140 (for NGC 5055). 
Comparing the disk-averaged extinction from this 
model to the radio-derived extinction, we find that 
for several galaxies, especially NGC 4321 and 5055, 
this method severely overestimates the global ex¬ 
tinction, by up to 4 mag. 

3. A simple alternative to the foreground dust screen 
model is a model in which the stars, gas, and dust 
are uniformly mixed. This model is motivated by 
the expectation that Ha traces stars which have re¬ 
cently formed from molecular gas, so the Ha and 
dust are likely to be physically co-extensive. For 
this case we have (e.g., Lequeux et al. 1981): 

Ay = 2.5 log - ——(12) 

where r is the extinction optical depth, again de¬ 
rived from the standard Galactic extinction law, 
t = Ay/1.086 = -/Vh[cui -2 ]/(2.17 x 10 21 ). In the 
limit of small optical depth, this reproduces the fore¬ 
ground screen model, Ay = 2.5 loge T / 2 , whereas 
for large r, Ay increases oc logr rather than oc r. 
This leads to much smaller extinction corrections at 
high gas densities, and closer agreement with the 
radio-derived extinctions. However, it appears to 
systematically underestimate the extinction for the 
galaxies with large Ay (e.g., NGC 5033). 

4. Finally, we consider a “hybrid” model in which the 
H 2 is mixed with the stars (Eq. 12) but the H I 
acts as a foreground screen, motivated by the larger 
scale height of the H I layer. When globally av¬ 
eraged, this leads to reasonable agreement (within 
~0.5 mag) with the radio-inferred extinctions, as 
shown in Figure 5(b). 

For the analyses in this paper we will employ both meth¬ 
ods (1) and (4). Method (1), a uniform correction across 
the galaxy, has the advantage that it is based directly on 
measurements (of Ha and radio fluxes) and agrees with 
previous studies that have found at most a weak trend in 
H II region extinction with galactocentric radius in M51 
and NGC 6946 (e.g., van der Hulst et al. 1988; Belley & 
Roy 1992; Hyman et al. 2000). Method (4), which pre¬ 
supposes a correlation between gas column density and 
extinction, agrees roughly with the radio fluxes but is not 
based on actual extinction measurements. If incorrect, 
it may lead to a spurious relationship between the gas 
density and star formation rate—hence the importance of 
applying both methods. The results of applying meth¬ 
ods (1) and (4) to the radial SFR profiles are shown by 
the short-dashed and long-dashed lines respectively in the 
lower panels of Fig. 3. 

A recent study by Quillen & Yukita (2001), based on a 
comparison of Ha with HST Paa imaging of a number of 
galaxies, including the CO-rich spirals NGC 2903, M51, 
and NGC 6946, lends some support to method (4). They 
find that the foreground screen model accurately predicts 
Ah q at l°w gas surface densities but overpredicts it at high 
densities. This is exactly the type of behavior provided by 
our hybrid model, since the Iow-Er regions are dominated 


by H I whereas high-Eu regions are dominated by H 2 . 
Moreover, the typical extinctions in the central regions of 
their galaxies are Ay ~ 2-4, consistent with the radial 
extinction profiles derived using method (4). It is worth 
noting that NIR recombination-line studies seldom reveal 
bright H II regions suffering extinctions of Ay > 4, even 
in circumnuclear starbursts (e.g., Kotilainen et al. 2000) 
or in situations where much larger extinctions might be 
expected based on the gas surface density (Krabbe et al. 
1994). 

Our extinction corrections may still underestimate the 
SFR if a substantial fraction of Lyman continuum pho¬ 
tons escape from galaxies or are absorbed by dust within 
the H II regions themselves, since such photons will not 
lead to observable Ha or radio continuum emission. The 
first effect, escape from the galaxy as a whole, is likely to 
be small -for example, Dove & Shull (1994) estimate an 
escape fraction for the Milky Way of ~14% (but see Beck¬ 
man et al. 2000 for a contrasting view). While this frac¬ 
tion may vary with radius, Ferguson et al. (1996) find little 
variation in the diffuse Ha fraction with radius, a quantity 
which they argue measures leakage of Lyman continuum 
photons. As for the second effect, dust absorption, Mathis 
(1986) estimates a typical optical depth to Lyman contin¬ 
uum photons in H II regions of r ~ 1.4, in which case we 
have underestimated the SFR by up to a factor of ~4. An 
analysis by McKee & Williams (1997), however, indicates 
that such absorption is less significant in H II region en¬ 
velopes than within the Stromgren radius; globally they 
find only ~25% of ionizing photons are absorbed by dust, 
yielding a correction factor of 1.37. Given the present lack 
of consensus, we apply no further corrections but note that 
our SFR estimates may be somewhat low. 

4. RESULTS 

4.1. The Star Formation Law 

Using the radial profiles of CO, H I, and total gas den¬ 
sity derived in §3.2, we can compare the gas surface density 
with the observed SFR surface density (derived from the 
Ha profile) on a ring-by-ring basis. Such a comparison 
involves averaging over greater areas for outer rings as op¬ 
posed to inner rings, and thus both random and systematic 
errors will vary with radius. In particular, the lack of inde¬ 
pendent data points may be a problem in the inner rings, 
whereas errors in flat-fielding or masking may dominate 
for outer rings. To minimize the impact of our masking 
method, we have only taken the radial profiles where they 
exceed the 3cr es t level defined in Equation 6. We note that 
despite their limitations, radial profiles are fairly insen¬ 
sitive to angular resolution and permit one to effectively 
average over time, which can be useful for revealing over¬ 
all trends. A point-by-point comparison of gas and star 
formation would be quite sensitive to the map resolution 
and to evolutionary effects, since on small enough scales, 
H II regions will disrupt their natal clouds. In a future 
paper, we will consider additional techniques for analyz¬ 
ing the two-dimensional intensity maps for both E gas and 
Esfr provided by BIMA SONG. 

4.1.1. The SFR and Molecular vs. Atomic Gas 

We first consider the relationships between Esfr and 
the atomic and molecular surface densities, Ehi and Er 2 , 
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Fig. 6.— SFR plotted against H 2 gas density (including associated helium), with both quantities azimuthally averaged, for the uniform 
extinction model. Dashed lines represent lines of constant star formation efficiency, with the corresponding gas depletion time labeled (as¬ 
suming H 2 is the only fuel for star formation). Note that we have taken advantage of the higher resolution (~6 ;/ ) of the CO data in making 
the comparison. 
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Fig. 7.— SFR plotted against H I gas density (including associated helium), with both quantities azimuthally averaged, for the uniform 
extinction model. Open symbols represent points where the molecular fraction is <0.5. Dashed lines represent lines of constant star formation 
efficiency, with the corresponding gas depletion time labeled (assuming H I is the only fuel for star formation). 
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separately. One advantage of such comparisons is that 
they avoid ambiguities resulting from the extrapolation of 
the CO profile; moreover, the comparison with £h 2 can be 
made at the higher resolution of the BIMA data. As de¬ 
scribed in §3.3, we consider two ways of correcting the Ha 
data for extinction: applying a uniform correction across 
the disk, or a radially varying correction that depends on 
the average gas column density. While the latter method 
may seem more physically motivated, evidence for radial 
extinction gradients is inconclusive at present. The results 
for the uniform extinction model are shown in Figures 6 
and 7, with up to two galaxies shown per panel. 

Most galaxies exhibit an excellent correlation between 
azimuthally averaged Ha and CO emission that can be 
parametrized using a Schmidt (power) law: 

_ dlogSsFR /"i q\ 

^-mol ,, v i (13) 

d log 5_]mol 

For the An-dependent extinction corrections, the inferred 
power-law indices range from n mo i « 1 for NGC 4414 to 
rimoi ~ 2 for NGC 5033 (see Table 5). The two excep¬ 
tions are NGC 4321, where the SFR is suppressed in the 
vicinity of a bar (i?=15"-55", where log£ mo i « 1.4), and 
NGC 4736, where the SFR is enhanced near a resonance 
ring (i?=30"-55", where log£ mo i « 1.6). These exam¬ 
ples suggest that the molecular content alone does not 
always control the SFR.—dynamical effects can also be 
important (see e.g., Wong & Blitz 2000; Meier & Turner 
2001). Even in NGC 4321 and 4736, however, there are re¬ 
gions where a power-law relationship is roughly obeyed. A 
weighted average for all seven galaxies gives n mo i=0.8 for 
the uniform extinction model and h mo i=1.4 for the vary¬ 
ing extinction model. For comparison, Rownd & Young 
(1999) made point-by-point CO and Ha comparisons for 
over 100 galaxies at 45" resolution, and found n mo i=1.2- 
1.4 with no extinction corrections applied. However, they 
note that n mo i=l (constant SFE for the molecular gas) 
is not excluded by their data. Assuming Ay does not 
vary strongly with radius, our study also finds reasonable 
agreement with a law in which the SFR per unit area is 
proportional to the molecular gas surface density. 

The relationship between Esfr and Ehi is more com¬ 
plex and appears inconsistent with a Schmidt law. In some 
cases, there is an abrupt decline in Ehi even when SgpR 
is fairly high (NGC 4736 and 5457, and to a lesser extent 
NGC 4414 and 4501). These galaxies have central H I de¬ 
pressions where significant star formation is still occurring. 
In other cases, the value of Ehi stays fairly constant over 
a wide range in Esfr, suggesting that a maximum value 
of Ehi has been reached. This maximum value, roughly 
10 Mq pc~ 2 or Ahi ~ 10 21 atoms cm -2 , could result from 
a tendency for H I to convert to H 2 at higher column den¬ 
sities due to self-shielding (Federman et al. 1979; Shaya 
& Federman 1987). Alternatively, the optical depth in 
the 21-crn line could become substantial above this col¬ 
umn density, especially if a cold H I component dominates 
(e.g., Braun 1997). The opacity at the line center is given 
by (Dickey & Lockman 1990): 

r = 5.2x 1 0-“(2|L) , (14) 

which produces r=l at Ahi=10 21 cm~ 2 for T=50 K and 
Au=10 km s -1 . However, Fig. 3 suggests that the max¬ 
imum deprojected Ahi is not dependent on inclination 


as one would expect if it were an optical depth effect: 
NGC 5033 (cos *=0.37) and 5457 (cos *=0.93) both exhibit 
azimuthally averaged H I surface densities that plateau 
around ~6-8 M 0 pc -2 . 

Surprisingly, Kennicutt (1989, 1998a) has found that 
the disk-averaged SFR is more strongly correlated with 
the disk-averaged H I surface density than with the CO 
surface density. Similar conclusions were drawn by Dehar- 
veng et al. (1994), who traced the SFR using far-ultraviolet 
(FUV) data, and Boselli (1994), who used both FUV and 
Ha. Our study, based on fully-sampled CO maps, is not 
consistent with these findings. The lack of a correlation 
between Esfr and Ehi in our data may result from a bias 
in our sample towards molecule-rich galaxies, where the 
H I surface density tends to saturate at a value of ~10 
Mq pc -2 . Galaxies with substantially less molecular gas 
may show a wider range of Ehi • It is more difficult to un¬ 
derstand why the correlation between global CO flux and 
SFR found by Kennicutt and others should be so poor, 
but K98 notes that the correlation is improved when only 
metal-rich spirals (such as those in our sample) are consid¬ 
ered, an effect he attributes to variations in the CO-to-H 2 
conversion factor (X-factor). Furthermore, the inhomo¬ 
geneity of the CO distribution in most galaxies suggests 
that incorrect H 2 masses may be derived when the CO 
distribution is not completely mapped. 

4.1.2. The SFR and the Total Gas Surface Density 

Combining the CO and H I data (with the CO profile ex¬ 
trapolated to larger radii) enables us to compare the SFR 
surface density with the total gas surface density, which has 
been argued to be a more relevant quantity for star for¬ 
mation than Eh 2 alone (Kennicutt 1989; K98). Figure 8 
shows the radial profiles for all seven galaxies as plotted 
in the £ gas -EgFR plane for the two extinction models. Re¬ 
gardless of how extinction is treated, there is indeed a 
strong correlation between £ gas and Esfr within galax¬ 
ies, indicative of a Schmidt law (Eq. 1). As is clear from 
the results of §4.1.1, however, this correlation is driven en¬ 
tirely by the molecular component for these galaxies, as 
Esfr shows virtually no correlation with Ehi- We note 
that the situation may be different in galaxies where the 
ISM is predominantly atomic, and further studies should 
be aimed at exploring the connection between star forma¬ 
tion and the H I component in such galaxies. 

Also shown in Fig. 8 as a dashed line is the composite 
Schmidt law derived by K98, 

S SF r = 0-16 f ^ gas 2 ) M 0 Gyr -1 pc -2 . (15) 

V1M 0 pc 

(Note that the numerical coefficient has been adjusted to 
include helium in E gas .) Kennicutt’s parameterization is 
based on averages over the star-forming disks of 61 nor¬ 
mal galaxies and the inner ~1 kpc disks of 36 starbursts. 
Two features are evident from Fig. 8. First, the use of 
an An-dependent extinction correction produces a steeper 
dependence of £ gas on Esfr than a uniform correction, as 
would be expected (e.g., Buat et al. 1989). Second, most 
of the observed data points lie to the right of Kennicutt’s 
relation. Note, however, that many of the points at larger 
E gas hr Fig- 8 occur at inner radii, which are weighted 
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Fig. 8. — SFR plotted against total gas density for rings within all seven galaxies, after applying (a) uniform extinction corrections; (b) 
IVe-dependent extinction corrections. Parallel dotted lines represent lines of constant SFR per unit gas mass, with the corresponding gas 
depletion time labeled. The heavy dashed line is the global Schmidt law derived by K98. 


less when taking a global average. Averaging our data 
across the entire observed disk region leads to much bet¬ 
ter agreement with Kennicutt’s results (Figure 9), even 
though our outer radius is generally smaller than the opti¬ 
cal radius adopted by K98. Indeed, we find that the offsets 
between our disk-averaged values and Kennicutt’s points 
for the same galaxies in Fig. 9 can largely be attributed to a 
change in the normalizing radius, which shifts points along 
a line of unit slope for the same total SFR and M gas . We 
conclude that within the uncertainty of ±0.2 dex adopted 
by K98, our results are consistent with his. 

Table 5 gives the average Schmidt law index, 


rflog Esfr 

d log Sgas 


(16) 


within each galaxy as derived from an unweighted least- 
squares fit to the curves in Fig. 8. An average of the indices 
for all seven galaxies, weighted by the inverse of their vari¬ 
ances, gives n = 1.1 ± 0.2 for the uniform extinction model 
and n = 1.7 ± 0.3 for the An-dependent extinction model. 
We examined the effect of changing the radial CO pro¬ 
file in the region where it had been extrapolated: in the 
extreme cases of a sharp truncation or complete flatten¬ 
ing of the CO profile, we obtained values of n = 0.8 and 
1.45 respectively for the uniform extinction model. How¬ 
ever, such extreme profiles seem unlikely, and we consider 
the la errors on n stated above to be realistic. Given 
the uncertainties, our derived indices are roughly consis¬ 
tent with the Schmidt law among galaxies (n ~ 1.4) found 
by K98. Note that there would not necessarily have been 
a correspondence between the azimuthally averaged and 
global Schmidt laws if star formation in galaxies depended 
on some quantity other than £ gas , but for which a disk- 
averaged (£ gas ) was a convenient proxy. Thus, the validity 
of the Schmidt law on both local and global scales proba¬ 
bly reflects an underlying physical link between £ gas and 
£sfr- 
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Fig. 9.— Disk-averaged SFRs plotted against total gas density. 
Large plot symbols are averages derived from radial profiles given 
in this paper, with radially varying extinction corrections applied. 
When connected by a line to a smaller symbol, the small symbol 
represents the corresponding point for the same galaxy as given by 
K98. The small dots are the other disk averages calculated by K98, 
the heavy dashed line is his global Schmidt law, and the error bars 
in the upper left represent his estimate of the uncertainties. 
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Table 5 

Observed Schmidt Law Indices for Total Gas and H2 Only 


Galaxy 

total 

Ay const. 

gas 

A v (N h ) 

molecular gas 

Ay const. Ay(Nn) 

NGC 4321 

1.13T.12 

1.74T.14 

1.41T.16 

2.08T.16 

NGC 4414 

1.15T.07 

1.71T.05 

0.32T.06 

0.98T.06 

NGC 4501 

0.95T.04 

1.60T.03 

0.90T.10 

1.63T.10 

NGC 4736 

0.82T.09 

1.23T.09 

0.04T.24 

0.74T.24 

NGC 5033 

1.38T.04 

2.06T.04 

1.27T.13 

2.01T.13 

NGC 5055 

1.14T.07 

1.79T.06 

0.55T.08 

1.30T.08 

NGC 5457 

1.18T.17 

1.56T.17 

0.88T.03 

1.45T.04 

Weighted Avg. 

1 .12T.22 

1.75T.25 

0.78T.34 

1.36T.35 


4.1.3. Gas Depletion vs. Orbital Timescales 

We have also considered an alternative description of 
the star formation law given by 

SsFR OC £ gas fl , (17) 

where H is the orbital frequency (e.g., Larson 1988; Silk 
1997). This law can be expressed as 

Tori) — CTg as , ( 18 ) 

where T gas = £ gas /£sFR the timescale for gas deple¬ 
tion by star formation, r or b is the orbital period and e is 
an efficiency factor that represents the fraction of avail¬ 
able gas that is consumed in each orbit. K98 finds that 
Equation 18, with e ~ 0.1, provides an equally good fit to 
the disk-averaged data in comparison to the Schmidt law. 
However, this prescription can only be properly tested us¬ 
ing azimuthally averaged data within galaxies, since a disk- 
averaged (t 0 rb) may simply be a reflection of the galaxy’s 
dynamical mass. 


Figure 10 shows the relation between r or b and r gas for 
the two extinction models we are considering. r or b was 
determined from rotation curves fitted to the CO and H I 
velocity fields, which had been derived from the datacubes 
by an automated Gaussian fitting routine (details will be 
given in Paper II). Only points where R > 30" have been 
plotted, since inside of that region the adopted rotation 
curve has a higher spatial resolution than the total gas 
profile. Applying a uniform extinction, we find that al¬ 
though a roughly linear trend among galaxies is apparent, 
as found by K98, no strong linear relation exists within 
galaxies. This is not surprising, since T gas tends to be 
constant with radius when uniform extinction corrections 
were applied, whereas Eq. 18 predicts r gas oc R for a flat 
rotation curve. The correlation improves somewhat for 
the varying extinction model [Fig. 10(b)], and indeed the 
scatter in log[EsFR/(S g as^)] is comparable to the scatter 
in log(£sFR/5!gat) f° r this model. Thus a star formation 
law where T gas oc r or b appears to be as valid a prescription 


Uniform A v N H -dependent A v 




Fig. 10.— Orbital period r or b = 2tr/fl plotted against gas depletion time r gas for annuli within each of the seven galaxies, after applying (a) 
uniform extinction corrections; (b) Ain-dependent extinction corrections. Parallel dotted lines denote points where 1% or 10% of the available 
gas is consumed by star formation during each orbit. 
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for the SFR as the Schmidt law, but only if the extinc¬ 
tion increases towards the center to allow for a significant 
gradient in T gas . To truly distinguish between the two pre¬ 
scriptions will require better constraints on the extinction, 
or the use of star formation tracers impervious to extinc¬ 
tion. 


4.2. Star Formation Thresholds 


Is star formation suppressed below a threshold or “crit¬ 
ical” gas density? The observation that H I in spiral 
galaxies often extends well beyond the optical disk has led 
to widespread support for the notion of a star formation 
threshold that is constant or depends on local conditions. 
As proposed by Quirk (1972), such a threshold may be re¬ 
lated to the parameter Q g for gravitational instability in 
the gas layer (Safranov 1960; Toomre 1964): 


Q g = , where E crit = ^77 ■ (19) 

Here n is the epicyclic frequency and c g is the gas veloc¬ 
ity dispersion. For a thin gaseous disk, values of Q g > 1 
correspond to stability and Q g < 1 to instability. Kenni- 
cutt (1989) found that in a sample of 15 mostly late-type 
galaxies, massive star formation was limited to regions of 
the disk where Q g < 1.6, suggesting that gas beyond some 
threshold radius is no longer susceptible to gravitational 
instability. Within the star-forming part of the disk, he 
also found that Q g did not vary much as function of ra¬ 
dius, suggesting that star formation is self-regulated so 
that Q g ~ 1 . 

We constructed radial profiles for E cr it using the 
adopted rotation curves presented in Paper II. The 
epicyclic frequency k is given by 


K 


2 


2 v c (dvc_ Vc\ 
~R \dR + ~Rj 


( 20 ) 


Note that for a flat rotation curve, n 2 = 2S2 2 . The deriva¬ 
tive of the rotation curve was determined by simple two- 
point interpolation. Following Kennicutt (1989), we as¬ 
sumed a value of 6 km s -1 for the velocity dispersion c g , 
which is consistent with detailed observations of H I in the 
outer regions of NGC 1058 (van der Kruit & Shostak 1984; 
Dickey et al. 1990) and of CO in NGC 628 (Combes & 
Becquaert 1997). The resulting profiles of E cr it are shown 
in Figures 11 12, along with profiles for E gas and Esfr. 
Note that E cr ; t oc i ? -1 for most galaxies, a reflection of the 
relatively flat rotation curves (D oc R~ 1 ). In a few cases 
(e.g., NGC 4736 and 5457), large spikes in E cr jt are seen; 
these are likely due to errors in determining dv c /dR via 
interpolation. 

At first glance, it is impressive how well the relation 
E gas « E cr jt holds for some of the galaxies (NGC 4501, 
5033, and 5055). On the other hand, there are two galax¬ 
ies (NGC 4321 and 4414) where E gas > E cr ; t over a sig¬ 
nificant range of radius, while in NGC 4736 and 5457, 
E gas < E cr ; t across most of the disk. The discrepancies are 
generally within a factor of 2 , which is certainly within the 
plausible range of error for either E gas (where a constant 
X-factor has been used and the CO profile extrapolated) 
or E cr j t (which depends on the assumed velocity disper¬ 
sion, inclination, and distance). Note, however, that the 
greatest uncertainties are probably associated with the X- 
factor and c g , and given our choices for these parameters 


we have probably overestimated E gas and underestimated 
Edit- Thus, the large inferred values for Q g = E cr ;t/E gas 
are almost certainly significant for NGC 4736 and 5457, 
yet star formation is hardly suppressed in these galaxies. 
(Surprisingly, these two galaxies are among those used by 
Kennicutt to argue for a Q threshold.) 

Indeed, we find no clear evidence that E cr it is relevant to 
Esfr hr the inner disk, since in most cases the E cr ; t pro¬ 
file is relatively featureless, a reflection of the relatively 
featureless rotation curve. For instance, in NGC 4321 it 
appears to be the case that star formation is enhanced 
when E gas > E cr i t and vice versa. Yet this may simply 
be due to the SFR profile roughly tracking the gas den¬ 
sity profile (see also Fig. 3). Similarly, the starburst ring 
in NGC 4736 is associated with an increase in E gas , not 
a dip in E cr ; t . Thus one should carefully examine claims 
that star formation occurs in a certain region because the 
gas density exceeds the critical density. The fact that the 
gas surface density is higher is probably relevant, but the 
critical density itself may not be. 

While the Q model was primarily intended to explain 
the outer disk cutoff of star formation (Kennicutt 1989), 
our data do not provide strong constraints in this regime. 
Deep, wide-held Ha imaging is required, ideally with spa¬ 
tially resolved measurements of the gas velocity disper¬ 
sion. Such a study has been undertaken by Ferguson et al. 
(1998), who found that a threshold radius could not be 
defined for NGC 6946 when using measured values for 
c g from Kamphuis (1993), although qualitative agreement 
with the Q model could be achieved if c g were assumed to 
be constant with radius. Since c g is very difficult to mea¬ 
sure in galaxies that are not face-on, it is unclear whether 
this result indicates a failure of the Q model or reflects 
inaccuracies in the measurement of c g . 

4.3. A Relationship Between Atomic and Molecular Gas 

The results of §4.1.1, which showed that the SFR traces 
the radial CO distribution and not the H I, imply that 
the SFR on large scales will be strongly influenced by 
the fraction of interstellar matter in molecular form. In 
a steady-state configuration, this fraction is determined 
by a balance between molecular cloud formation and de¬ 
struction processes. The timescales of these processes and 
the amount of material involved can only be determined 
empirically when there are organized features, such as spi¬ 
ral arms, with which one can trace the H I/H 2 ratio az- 
imuthally. In the inner disk of M51, for example, virtually 
all of the H I is believed to result from dissociation of H 2 in 
the spiral arms (Vogel et al. 1988; Tilanus & Allen 1989; 
Rand et al. 1992). Conversely, Heyer & Terebey (1998) 
find that the gas entering the Perseus arm in the outer 
Galaxy is primarily atomic; in this case the H I must be 
a precursor to the H 2 in the arm. We will defer discus¬ 
sion of time-dependent processes to a future paper that 
includes azimuthal comparisons, and will focus instead on 
what determines the time-averaged H I/H 2 ratio. Thus we 
ignore the question of whether H I is a product of H 2 or 
vice versa, since in a time-averaged sense the distinction 
has little meaning. 

In comparing the azimuthally averaged values of Ehi 
and Eh 2 , we are restricted to a fairly narrow region in the 
disk: R > 20", in order to minimize the effect of beam 
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NGC 4321 NGC 4414 




NGC 4501 NGC 4736 




Fig. 11.— Comparison of total gas density with the critical density for gravitational instability for NGC 4321, 4414, 4501, and 4736. The 
Ha profile, uncorrected for extinction, is also shown for comparison (dotted line). The error bar at the lower left represents a factor of 2 
change, and the heavy dashed line represents an R~ 1 profile, as expected for E cr it with a flat rotation curve. 
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NGC 5033 NGC 5055 




NGC 5457 NGC 5457 




Fig. 12.— Same as Figure 11, for NGC 5033, 5055, and 5457. For NGC 5457 the horizontal axis is also shown on a linear scale (lower 
right ). 
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smearing on the profile, and R < 80", the extent of our 
CO observations. Furthermore, we have rejected points 
where E H i < 1 M 0 pc -2 , because the FI I profile is gen¬ 
erally uncertain at these radii (this was only an issue for 
R < 36" in NGC 5457). In practice, this means that we 
only probe a region where the molecular fraction j mo \ falls 
from ^0.9-1 to ~0.6-0.7, emphasizing the need for high 
sensitivity, large-field CO data to fully analyze the transi¬ 
tion region. Figure 13 compares the azimuthally averaged 
H I and H 2 surface densities for all seven galaxies over the 
region satisfying the above criteria. The two types of be¬ 
havior seen in the radial gas profiles in Fig. 3 are reflected 
in whether Ehi is constant with increasing Eh 2 or begins 
to decrease at some point. In NGC 4414, 5033, and 5055, 
the H I profile remains relatively flat across the inner disk, 
while in NGC 4321, 4501, 4736, and 5457, the H I profile 
declines toward the center. 



Fig. 13.— Relation between molecular and atomic gas surface 
density for all seven galaxies in the region where observations over¬ 
lap (Rmin=20 // , Rmax ~ 80") and where Shi > 1 M® pc -2 . 


As shown in Figure 14, the ratio of Ehi to Eh 2 appears 
to depend on radius, increasing as roughly 

= KR 15 (21) 

£h 2 

over a factor of ~2 (0.3 dex) in R. The fact that the curves 
in Fig. 14 are parallel but are spread out over nearly an 
order of magnitude in R indicates that there is a second 
parameter K that is roughly constant within a galaxy but 
varies from galaxy to galaxy. If K is taken to be (i? 2 s) -1 ' 5 , 
where R 25 is the optical radius given in Table 1, then the 
residual variation in logi?, although still substantial, is 
reduced to ~0.5 dex. The most striking aspect of the re¬ 
lation given by Equation (21) is that it is independent of 
whether the H I profile is flat or declining towards the 
center: note especially how it “straightens out” the lines 
for NGC 4321, 4736, and 5457 when compared to Fig. 13. 
Thus, the ratio of Ehi to Eh 2 appears to depend on galac- 
tocentric radius, normalized by R 25 , and little else. In 
galaxies where the inner H I profile is flat, Eh 2 declines 


as roughly i? -1 - 5 , whereas the galaxies with H I decreas¬ 
ing towards the center have a correspondingly flatter H 2 
profile across the transition region. We discuss a possible 
interpretation of these trends in §5.4. 



Fig. 14.— Ratio of H I to H 2 surface density plotted against 
radius for all seven galaxies. Two power-law parameterizations are 
shown, one where the ratio scales with R and another where the ra¬ 
tio scales with R 1,5 . The data are the same as shown in Figure 13. 


5. DISCUSSION 

5.1. Does CO Trace Molecular Mass? 

We have interpreted the linearity between azimuthally 
averaged CO and Ha emission (§4.1.1) as indicative of a 
constant star formation efficiency, but it could conceivably 
be a result of CO excitation effects. That is, if energy input 
from young stars is responsible for exciting the CO emis¬ 
sion, then CO may trace star formation directly rather 
than the H 2 mass (e.g., Dopita & Ryder 1994). Evidence 
in favor of this viewpoint comes from observations of the 
[C II] 158 pin line, a tracer of photodissociation regions, 
which is found to scale linearly with the CO(l-O) flux in 
infrared-luminous galaxies (Crawford et al. 1985). How¬ 
ever, relatively stronger CO emission (compared to [C II]) 
is found in spirals with more typical SFRs (Stacey et al. 
1991), and further studies with higher-resolution instru¬ 
ments such as SOFIA will be essential. 

In any case, there are at least three difficulties with the 
suggestion that we have significantly overestimated the H 2 
mass in regions of active star formation. First, the gas de¬ 
pletion times are already quite short in these regions, ~1 -3 
Gyr, and become even shorter if there is less H 2 ; it seems 
unlikely that we would happen to witness the final stages 
of star formation in these regions. Second, observations 
of cold dust in nearby galaxies with the SCUBA submil¬ 
limeter camera show a strong correlation between CO and 
850 pm emission (e.g., Alton et al. 2001), indicating that 
variations in the X-factor with radius cannot significantly 
exceed variations in the gas-to-dust ratio. Finally, the spa- 
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tial correspondence between CO and young massive stars 
(which would dominate the radiative excitation) does not 
appear to be close when observations of sufficient angular 
resolution are performed. Note for instance the different 
azimuthal symmetries of CO and Ha in NGC 4736 (Wong 
& Blitz 2000), or the spatial offset between the brightest 
CO emission and the most populous young clusters in the 
LMC (Cohen et al. 1988). 


5.2. Interpretation of the Schmidt Law 


Previous evidence for a power-law relationship between 
Ega S and Esfr, with n « 1.5, has been discussed by Ken- 
nicutt (1997) and K98. As shown in §4.1, this relation 
appears to hold within galaxies as well. The most straight¬ 
forward interpretation of such a Schmidt law is that the 
SFR on large scales is controlled by the self-gravity of the 
gas. As a result, the rate of star formation is proportional 
to the gas mass divided by the timescale for gravitational 
collapse: 

gas 1.5 / 0 o\ 

dSFR — jg j-0.5 ^ Pgas • (22) 


The problem with this explanation, however, is that the 
efficiency factor e* is observed to be much less than 1. 
In other words, the gas depletion time on large scales 
(jgas ~ 10 9 yr), while often much less than a Hubble time, 
is still significantly longer than Tg = (Gp gas )~ 0 ' 5 « 10'- 
10 s yr for «h ~ 1- 100 cm -3 . This leads to two competing 
interpretations: (1) the star formation timescale really is 
Tg, but only a small fraction e* of the gas participates in 
star formation (Elmegreen 2000), or (2) the star formation 
timescale is Tg multiplied by a factor of e/ 1 , due to the 
inhibiting effects of stellar feedback (i.e., winds, radiation, 
and supernovae) or magnetic support on star formation 
(Larson 1988). Neither description, however, offers a clear 
prediction for the value of e* or explains why it should 
be constant. Replacing rg with the orbital timescale T or b> 
as discussed in §4.1.3, yields a similar problem: why is a 
small but fixed fraction of the gas consumed per orbit? 

Given that the theoretical basis of the Schmidt law re¬ 
mains unclear, it is worthwhile to consider an alternative 
explanation, motivated by our finding in §4.1.1 that Esfr 
scales roughly linearly with £h 2 —a result that is also con¬ 
sistent with previous extragalactic (Rownd & Young 1999) 
and Galactic (McKee 1989) studies. In this picture, stars 
form only in molecular gas, and the timescale for star for¬ 
mation is controlled by stellar feedback, magnetic diffu¬ 
sion, or other processes which do not vary strongly with 
location in a galaxy, so that Esfr oc Eh 2 (n mo i=l). As 
suggested by Tacconi & Young (1986) and Wang (1990), 
the Schmidt law between Esfr and the total gas surface 
density E gas follows from an amalgamation of two “laws,” 
the first relation governing the formation of H 2 from H I, 
and the second (linear) relation governing the formation of 
stars from H 2 . The composite Schmidt law index is given 

by 


d In Esfr 
dln£ 


gas 


/ din Esfr A ( din E H 
V din£ h 2 ) 
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For n mo i « 1 and f mo \ increasing with £ gas (since 
/moi increases with interstellar pressure which scales with 
Egas—see §5.4 below), one naturally arrives at n > 1. 
Physically this also makes more sense than a star forma¬ 
tion law based on the total gas surface density, since when 
one looks in the third dimension, one finds that the ver¬ 
tical distribution of star formation corresponds to that of 
H 2 , not H I. The sharp increase in the Schmidt law index 
n seen at low E gas , which Kennicutt (1989) attributes to 
a Q-dependent threshold for star formation, could instead 
be due to a sudden decrease in / mo 1 at low E gas (or low 
p gas , if the volume density is more relevant due to flar¬ 
ing of the HI layer, Ferguson et al. 1998), although the 
two thresholds may turn out to be equivalent if molecu¬ 
lar cloud complexes are formed by large-scale gravitational 
instabilities. 

This “two-step” model is not without its difficulties—it 
still does not predict the SFE in molecular clouds quan¬ 
titatively, nor does it explain the enhanced SFE’s found 
in merging galaxies (Young et al. 1996) or starburst rings 
(Wong & Blitz 2000), which could be related to cloud colli¬ 
sions (e.g., Tan 2000). However, it does provide an attrac¬ 
tive explanation for the relatively constant star formation 
histories of disks (Kennicutt et al. 1994), namely that the 
requirement of converting H I into H 2 could limit the SFR 
at early times, despite a greater supply of gas. Rana & 
Wilkinson (1986) have presented such a model, in which 
the H 2 formation rate at early epochs is limited by the 
lack of metals, since metals serve as important coolants 
in the ISM and constitute the dust grains on which most 
of the H 2 forms. Even today, the confinement of H 2 to a 
thinner disk than the H I (see e.g., studies of NGC 891 
by Scoville et al. 1993; Swaters et al. 1997) suggests that 
the formation of H 2 is inhibited outside of the midplane, 
perhaps because of insufficient pressure or dust abundance 
(see also the discussion in §5.4 below). 


5.3. Is Q Relevant to Star Formation? 

We have found that the Q g parameter, defined using a 
constant velocity dispersion of 6 km s , appears to be 
remarkably close to 1 in some galaxies, but is significantly 
greater than 1 in NGC 4736 and 5457, which are still ac¬ 
tively forming stars. Our observations thus lend support 
to earlier studies which have failed to show a clear relation¬ 
ship between Q g and the presence of massive star forma¬ 
tion (Thornley & Wilson 1995; Meurer et al. 1996; Fergu¬ 
son et al. 1998; Hunter et al. 1998; Pisano et al. 2000). This 
is not altogether surprising; Elmegreen (1999) raises the 
question of “why Q < 1 has anything at all to do with star 
formation in a supersonically turbulent medium filled with 
supernovae, magnetic Helds, multiple gas phases and other 
seemingly important details.” Yet the problem remains of 
explaining why Q g ~ 1 is often observed in galaxy disks, a 
fact that has long been taken to indicate that star forma¬ 
tion on large scales is regulated by gravitational instability 
(Quirk 1972; Kennicutt 1989; Silk 1997). A possible clue 
is that Q g is significantly greater than 1 in the two galax¬ 
ies with the lowest gas fractions (NGC 4736 and 5457), as 
inferred from a comparison of E gas with the stellar surface 



Gas and Star Formation in Spiral Galaxies 


25 


density £* derived by assuming a maximum disk (see Pa¬ 
per II). This suggests that Q g is actually a measure of the 
gas fraction in the disk. 

A simplified analysis suggests that this is indeed the 
case. For a Mestel (1963) disk, with £ tot a i? _1 , a flat 
rotation curve is obtained such that 


v\ — 2-7rG£tot-R > 


v c Tl 

= 27rG ' 


(27) 

(28) 


[See Binney & Tremaine (1987), p. 76.] But for a flat 
rotation curve, the critical density is given by 
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crit — 


\/2 flCg 

7 tG 


Defining the gas fraction as p = £ gas /£tot, we have 
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Substituting typical values for spiral galaxies (c g =6 
km s _1 , i> c =200 km s _1 , and /i=0.1) we find that Q=0.84. 
Thus, the fact that Q g is of order unity in many spiral 
galaxies reflects the fact that the present-day gas fraction is 
roughly 0.1. Note that we have implicitly assumed a max¬ 
imum disk model with negligible bulge or halo, but this 
analysis will be roughly appropriate if Etot is considered 
to be whatever mass distribution dominates the rotation 
curve. For instance, dwarf galaxies exhibit large Q g values 
(e.g., Meurer et al. 1996; van Zee et al. 1997) because they 
rotate slowly and have gas surface densities that are small 
compared to the dark matter surface density. We therefore 
propose that Q g ~ 1 is primarily a reflection of observa¬ 
tional selection favoring luminous, gas-rich galaxies for this 
kind of analysis; galaxies with smaller gas fractions should 
show Q g > 1. We note that this explains the larger Q g 
values found in early-type galaxies by Martin & Kennicutt 
(2001). Still, the lack of galaxies with Q 9 < 1 suggests 
that strongly unstable gas disks are precluded. 

That the Q g parameter is most sensitive to the gas frac¬ 
tion and thus cannot be the same for all galaxies was rec¬ 
ognized many years ago by Larson (1988). He suggested 
that an increase in c g would partially offset the larger £ gas 
at early times, maintaining Q g ~ 1. It is unlikely, how¬ 
ever, that the velocity dispersion of the gas, when averaged 
over large scales, can continue to decrease much below a 
value of ~6 km s -1 , since measurements in the face-on 
galaxy NGC 1058 indicate that it does not fall below this 
value even well beyond the optical radius (van cler Kruit 
& Shostak 1984; Dickey et al. 1990). Sellwood & Balbus 
(1999) attribute this minimum velocity dispersion to mag- 
netohyclrodynamic turbulence that is ultimately fed by the 
galaxy’s differential rotation. Based on these arguments, 
a secular increase in Q g seems inevitable as gas consump¬ 
tion continues, unless new gas is constantly accreted to 
maintain p at a constant level. 

In spite of the apparent difficulties in using Q g as a 
predictor of star formation activity, it is still remarkable 
how well Q g ss 1 is satisfied in certain galaxies, such as 
NGC 5033 and 5055 (Fig. 12). A possible explanation is 
that gravitational instability does regulate star formation, 


but the appropriate Q parameter must take into account 
both the gas and stellar disks. Elmegreen (1995) and Jog 
(1996), extending the work of Jog & Solomon (1984), have 
shown that the gravitational interaction of these two com¬ 
ponents leads to greater instability in the combined system 
than in either component individually. They derive an “ef¬ 
fective” Q parameter Q e ff which is « Q g for gas-dominated 
disks, but can be close to Q* when the gas fraction is small. 
For Q* to be regulated at a value of ~1 requires only that 
the stellar velocity dispersion c* increase in proportion to 
the stellar surface density, and is in good agreement with 
observations if one assumes a constant M/L ratio in the 
disk (Bottema 1993). Thus, even galaxies with low gas 
fractions may be susceptible to mild gravitational insta¬ 
bilities as a result of Q* ~ 1, and gravitational instability 
may be able to account for the presence of continued star 
formation in these galaxies after all. 

As noted by Skillman (1996), an alternative to a de¬ 
based star formation threshold is a constant column den¬ 
sity threshold of IVh ~ 10 21 cm -2 associated with the 
formation of molecular clouds. In this picture, star forma¬ 
tion cuts off at a radius where there is insufficient self¬ 
shielding from UV radiation to allow molecular clouds 
to form. In practice, the column density threshold may 
vary with metallicity, pressure, and the interstellar radi¬ 
ation field, all of which affect the transition from H I to 
H 2 (Elmegreen 1993). A similar criterion, in which the 
minimum column density is defined by the requirement 
of sufficient pressure to allow a two-phase ISM, has been 
discussed by Elmegreen & Parravano (1994). Supporting 
this hypothesis, Braun (1997) has presented evidence that 
cold H I clouds are not observed outside the optical radius 
of galaxies; the inability to form such clouds may be re¬ 
sponsible for star formation thresholds. A comparison of 
sensitive CO, H I, and Ha data near the edges of galaxy 
disks would be needed to test these ideas. 

5.4. The Atomic-to-Molecular Gas Ratio 

In §4.3 we found that the ratio of £hi to £h 2 varies 
within galaxies as roughly R 1 5 in the region of overlap. 
However, the interpretation of this dependence is ambigu¬ 
ous, since many galaxy properties vary with radius. More¬ 
over, the limited range in radius over which the ratio is 
calculated may allow for a variety of functional forms, of 
which this relation is only the simplest. For purposes of 
discussion, we consider here two possible models for how 
£hi and £h 2 should be related, and compare them with 
our results. 

First, Wyse (1986) has suggested that the formation of 
GMCs is tied to collisions between atomic clouds as they 
pass through spiral arms. Assuming a negligible pattern 
speed (i.e., well inside the corotation radius), the rate of 
spiral arm crossings at radius R is oc D(i?), and one ex¬ 
pects 

£h 2 cx (£hi) 2 ^ • (32) 

This naturally introduces an R dependence into the ra¬ 
tio £hi/£h 2 , qualitatively similar to the R 15 dependence 
we derived above. However, Figure 15 shows that while 
the majority of the galaxies show some correlation be¬ 
tween £h 2 and (£hi) 2 H, the expected relation (given by 
the dashed line) is a much poorer description of the data 
than the simple R 15 relation shown in Fig. 14. 
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log (£g,Q) 

Fig. 15.— Surface density of molecular gas vs. the prediction from 
Wyse (1986), where Eh 2 ^hi 2 ^- The dashed line represents a 
slope of 1, in agreement with the prediction. 

Second, Elmegreen (1993) argues that the degree to 
which a cloud is molecular is determined in part by the 
ambient pressure P, since the molecule formation rate de¬ 
pends on the volume density n, which scales with P for 
both diffuse and self-gravitating clouds. The pressure at 
the midplane of a two-component disk in hydrostatic equi¬ 
librium is given by (Elmegreen 1989): 

P « |GE gas (s gas + , (33) 

where c g and c* are the velocity dispersions of the gaseous 
and stellar disk respectively. The observations of Bottema 
(1993) typically give c* ~ 60 km s _1 , so we assume a value 
for c s /c* of 0.1. (This value for c* is consistent with the 
velocity dispersion of an isothermal disk, c* = 
for E* ~ 300 M 0 pc -2 and zo ~ 1 kpc.) We have de¬ 
rived radial profiles of E* using K or /-band profiles with 
a mass-to-light ( M/L) ratio chosen to provide a maximal 
disk contribution to the rotation curve, as discussed fur¬ 
ther in Paper II. As shown in Figure 16(a), a relation of 
the form 

Ehi/E H2 oc P- 0 - 8 (34) 

appears to match the data reasonably well. For compari¬ 
son, Elmegreen (1993) predicts a very sensitive dependence 
of the molecular fraction / mo i on pressure: / mo i oc P 2 ' 2 j , 

where j is the volume emissivity of dissociating radiation. 
For j oc E* oc E gas , this leads to / mo i oc P 1 , which is 
clearly steeper than what is observed within our galaxies 
[Fig. 16(b)]. However, his parametrization is intended for 
use in the low / mo i regime, which is more appropriate for 
the Milky Way than for ^-dominated galaxies; indeed the 
match to the Galactic data of Dame (1993) is reasonably 
good [Fig. 16(b)]. 

Residual offsets between galaxies in Fig. 16(a) may re¬ 
sult from uncertainties in the adopted M/L ratio, or 
could be related to additional physical parameters, such 
as metallicity. Specifically, we find that galaxies with 
oxygen abundances significantly above solar in their in¬ 
ner (P < 200") disks (NGC 4321, 4501, 5055, and 5457) 


display smaller values of Ehi/Eh 2 - Higher metallicity is 
expected to favor the formation of H 2 , since dust atten¬ 
uates the dissociating UV radiation field and provides a 
surface on which to form H 2 . The importance of metallic¬ 
ity in the H I-H 2 balance is discussed further by Rana & 
Wilkinson (1986) and Honma et al. (1995). However, we 
caution that the use of CO to trace H 2 may also introduce 
a metallicity dependence in our derived value of Ehi/Eh 2 - 
Thus, while the range in f mo \ probed by our observa¬ 
tions is fairly small, the trend in the H I/H 2 ratio with 
radius suggests that pressure is playing the dominant role, 
with metallicity as a secondary factor. This is supported 
by observations indicating that the global Mhi/Mh 2 ratio 
increases from early-type to late-type galaxies (Young & 
Knezek 1989; Casoli et al. 1998), since early-type galaxies 
generally have higher mass surface densities (and metallic- 
ities). Moreover, the tendency for recent star formation to 
be correlated with the distribution of stellar light, noted 
for example by Ryder & Dopita (1994) and Hunter et al. 
(1998), can be explained if stars currently dominate the 
pressure which controls H 2 and star formation. Further 
observations spanning a wider range of / mo i would be use¬ 
ful to verify the importance of pressure in determining the 
H I/H 2 balance. 


6. SUMMARY AND CONCLUSIONS 

We have investigated the relationship between gas con¬ 
tent and star formation rate, commonly termed the star 
formation law, in seven nearby galaxies using CO, HI, and 
Ha images at resolutions of ~20" or higher. An impor¬ 
tant advance is the recent availability of CO data that 
combines the resolution of an interferometer with the sen¬ 
sitivity to large-scale structure afforded by a single-dish 
telescope. Although our sample is biased toward lumi¬ 
nous, molecule-rich galaxies, none of the galaxies display 
signs of strong interactions. 

We find that the correlation of the azimuthally averaged 
SFR surface density Esfr with Ehi is much poorer than 
with Eh 2 , contrary to the results of studies based on in¬ 
tegrated fluxes. Whereas a roughly linear relation exists 
between Esfr and Eh 2 , consistent with a constant star 
formation efficiency for the molecular gas, Ehi tends to 
reach a maximum value of ^10 M 0 pc -2 , or even decline 
in regions of high SFR. Thus the star-forming gas in these 
galaxies exists predominantly in molecular form, although 
the situation may be different in low-mass galaxies. 

Despite the poor correlation of Esfr with Ehi, the dom¬ 
inance of H 2 in these galaxies leads to a strong correlation 
between Esfr and total gas surface density E gas when az¬ 
imuthally averaged, in rough agreement with the global 
Schmidt law derived by K98. The index of the power law 
is given by n « 1.1 if the extinction at Ha is assumed to be 
uniform with radius, but steepens to n w 1.7 if the mean 
extinction is assumed to vary according to the gas surface 
density. A star formation law of the form Esfr oc E gas fi is 
also consistent with the data, but only if the .An-dependent 
extinction corrections are applied. We suggest that the ob¬ 
served Schmidt law results from an effective dependence 
of the molecular fraction on E gas by way of the interstellar 
pressure, coupled with a roughly linear relation between 
Esfr and E H2 - 
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Fig. 16. — (a) Ratio of H I to H 2 surface density plotted against the hydrostatic disk pressure at the midplane, as given by Elmegreen 
(1989). The straight dashed line represents the mean slope of —0.8. (b) Molecular fraction plotted against midplane pressure. The straight 
dashed line represents the slope of 1.7 predicted by Elmegreen (1993). Data for the Milky Way is based on Dame et al. (1993), adopted a 
stellar scalelength of 3 kpc (Sackett 1997) and E* = 35 Mq pc -2 at the solar circle (Gilmore et al. 1989). 


The gravitational instability parameter Q g is found to 
be ~ 1 in several galaxies but is clearly > 1 in the two 
galaxies with the smallest gas fractions, NGC 4736 and 
5457. We suggest that Q g is best considered a measure 
of the gas fraction, and that when discussing a range of 
galaxy types, the contribution of the stellar component 
to the disk instability should not be neglected. While a 
combined (gaseous and stellar) instability may ultimately 
provide a better description of star formation thresholds, 
alternative explanations that attribute such thresholds to 
the inability to form cold, dense clouds remain attractive 
as well. 

Finally, over the limited range in radius where our CO 
and H I data overlap, we find that the relation Ehi/Eh 2 ^ 
R 1 - 5 holds remarkably well within galaxies, although the 
proportionality constant varies among galaxies. We show 
that such a relation supports the view that the interstel¬ 
lar pressure plays the dominant role in determining the 
H I-H 2 balance. 

These results have demonstrated the importance of in¬ 
cluding the molecular component in studying the star for¬ 
mation law. Moreover, considering just the total gas den¬ 
sity rather than the H I and H 2 surface densities separately 
may obscure underlying physical processes that are essen¬ 
tial to star formation. In a future study, we will develop 
these ideas using the much larger CO database provided 
by the BIMA SONG project, supplemented by archival 
and recently obtained H I synthesis imaging. We also em¬ 
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phasize that significant advances will be made possible by 
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